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THEORETICAL MONOCHROMATIC-WAVE-INDUCED CURRENTS 


IN INTERMEDIATE WATER WITH VISCOSITY 
AND NONZERO MASS TRANSPORT 

Theodore A. Talay 
Langley Research Center 

SUMMARY 

This study reviews some wave -induced mass-transport current theories with both 
zero and nonzero net mass (or volume) transport of the water column. A relationship 
based on the Longuet-Higgens theory is derived for wave-induced, nonzero mass- 
transport currents in intermediate water depths for a viscous fluid. The relationship 
is in a form useful for experimental applications; therefore, some design criteria for 
experimental wave-tank tests are also presented. Sample parametric cases for typical 
wave-tank conditions and a typical ocean swell were assessed by using the relation in 
conjunction with an equation developed by Unluata and Mei for the maximum wave- 
induced volume transport. Calculations indicate that substantial changes in the wave- 
induced mass-transport current profiles may exist dependent upon the assumed net 
volume transport. A maximum volume transport, corresponding to an infinite channel 
or idealized ocean condition, produces the largest wave -induced mass-transport currents. 
These calculations suggest that wave-induced mass-transport currents may have consid-. 
erable effects on pollution and suspended-sediments transport as well as buoy drift, the 
surface and midlayer water- column currents caused by waves increasing with increasing 
net volume transports. Some of these effects are discussed. 

INTRODUCTION 

The waters over continental margins are assuming a greater importance to many 
nations as considerations of economy, pollution, and transportation pose new requirements 
on our general knowledge of these waters. One major factor in the study of these waters 
is the oceanic circulation pattern, both global and local. Often, the local ocean current 
pattern must be investigated to determine the drift of a particular oil spill, ocean spoils 
dump, buoy, or school of fish. This local ocean circulation pattern is complicated by 
many interacting physical processes such as the ocean density field, sea surface topog- 
raphy, wind pattern and sea state, bathymetry, tides, zones of convergence and divergence 
of oceanic waters, and proximities of coasts. 


A recent study (ref. 1) examined some of these conditions in a region of the mid- 
Atlantic continental shelf. One of the conclusions of this study was that currents caused 
by waves are often of the same order of magnitude as other currents in the region. 

Waves in passage are known to impose a nonclosed orbital motion on the fluid particles 
beneath the passing waves (ref. 2). Alternatively, this condition may be described as a 
closed orbital motion with a drift associated with the progression of the orbits. In 
shallow waters the orbital motion becomes flattened and resembles a back and forth 
"scrubbing" velocity which may lift sediment and wastes off the bottom and place them 
in suspension. The wave -induced drift may then act to transport these sediments or 
wastes. At or near the surface this wave-induced drift may influence the drift of sur- 
face pollutants and buoys, and thus represent a major component of the surface drift. A 
number of theories have been advanced to explain and predict this wave-induced drift, 
each theory limited by its own particular set of assumptions. Experimental wave-tank 
data in the form of drift velocity profiles have both confirmed the presence of a wave- 
induced drift and provided comparisons of the various theories. 

Stokes' theory (ref. 2), assuming an inviscid fluid, presents equations for the wave- 
induced mass-transport current profile for cases of zero and nonzero net wave -induced 
mass transport of the water column. From experimental data for the zero net mass- 
transport cases, Stokes' theory appears to work well in deep waters. However, this 
theory does not explain observed velocities, especially a strong forward flow at the 
bottom in intermediate waters which are typical of continental shelf areas. The theory 
of Longuet-Higgens (ref. 3) follows an Eulerian approach to the equations of motion 
and allows for a nonzero fluid viscosity. For cases of zero net wave-induced mass 
transport of the water column, the Longuet-Higgens equation for the wave-induced 
mass-transport current profile qualitatively explains many of the effects observed in 
wave-tank experiments, even when the equation is applied beyond its original limitations. 
Recently, Unluata and Mei (ref. 4) examined the wave-induced velocity profiles from a 
Lagrangian approach to the equations of motion. Equations for the wave -induced mass- 
transport current and the corresponding net volume transport of the water column are 
presented. These equations involve wave and depth conditions as well as a constant 
pressure gradient in the flow direction. For a case of zero net mass (or volume) trans- 
port, the current profile was shown to reduce to that obtained by Longuet-Higgens. A 
maximum net mass transport of the water cplumn was stated as that corresponding to a 
zero pressure gradient. 

Kinsman noted that the assumption of a zero net horizontal mass transport caused 
by waves is appropriate for wave experiments carried out in a closed wave tank (ref. 5). 
In the ocean, however, it is unlikely that the net horizontal mass transport caused by 
waves will be zero. 
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This report discusses wave-induced current profiles associated with nonzero wave- 
induced mass transports. It reviews some previously published mass-transport current 
theories with a zero and a nonzero net mass transport of the water column. The appendix 
contains a derivation following closely the Eulerian approach of Longuet-Higgens (ref. 3), 
which results in an expression for the wave-induced mass-transport current profiles for 
intermediate water with viscous flow allowing a net mass transport of the water column. 
This expression is shown to be the same as that proposed by Longuet-Higgens although 
he did not explicitly state it in the form presented here. For experiments conducted in 
wave tanks where only current profiles may be generally measured, the derived expres- 
sion is shown to be in a form directly applicable for comparing theoretical wave -induced 
currents with the measured data. Necessary wave-tank design criteria for such tests are 
discussed. In the absence of such experimental data, two parametric evaluations were 
conducted, one for wave and depth conditions from a typical wave -tank test and the other 
for a typical swell condition on the mid-Atlantic continental shelf. A theoretical expres- 
sion for the maximum nonzero mass transport, given by Unliiata and Mei in their 
Lagrangian analysis (ref. 4), was used in the analysis which follows to determine the 
upper limit for the net wave -induced mass (or volume) transport. A range of net mass 
(or volume) transport from zero to the theoretical maximum was used. 

SYMBOLS 




a 

C 


H 

h 

k 

L 

£ 


coefficients in equation for Uljj(z) (defined by eqs. (7)) 
wave amplitude, meters 

coefficients in equation for MT^^^ (defined by eqs. (A20a) and (A20b)) 

coefficient used in Stokes' mass-transport current (defined by eq. (3)), 
meters/second 

wave height, trough to crest, 2a, meters 
mean water-column depth, meters 
wave number, 27 t/L, meters"^ 
wavelength, meters 
wave -tank length, meters 
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0( ) order of 

second-order, time-averaged pressure gradient, newtons/meter^ 

R(z) function (defined by eqs. (A17) and (A18)) 

T wave period, seconds 

U(z) wave-induced mass-transport current, meters/second 

net horizontal wave-induced volume transport per unit length of crest, 
meter s^/second 

X spatial coordinate measured positive in direction of wave advance, meters 

z vertical coordinate measured positive downwards from mean water level. 



meters 


0-1 

coefficient of 

R(z), meters"^ 

“2 

coefficient of 

R(z), meters"^ 

"3 

coefficient of 

R(z), meters”^ 


coefficient of 

R(z) 

6 

boundary -laye 

r thickness, meters 


e small dimensionless quantity used in series expansion, 0(a/L) 

II dynamic viscosity of water, kilograms/meter -second 

u kinematic viscosity of water, meters^/second 

p density of water, kilograms/meter^ 

a wave frequency, 2tt/T, seconds"^ 

i// stream function for interior of fluid, meter-meter/second 
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Subscripts: 


B bottom 

DP in deep water 

EXP experiment 

LH Longuet-Higgens, under conditions of zero net mass transport 

MAX maximum 

MT conditions of nonzero net mass transport 

S surface 

UM Unluata and Mei, under conditions of zero net mass transport 

REVIEW OF SOME WAVE -INDUCED CURRENT THEORIES 

With the passage of a water wave, the fluid particles beneath the wave surface 
theoretically exhibit a nonclosed orbital motion in the vertical plane as shown in the 
following sketch: 




+x 


f +z 

This condition may be described, alternatively, as a closed orbital motion with a drift 
associated with the progression of the orbital motion. This steady, second-order drift 
velocity (also known as the mass-transport current) has been confirmed by several 
experimenters (ref. 3, p. 577; refs. 6 and 7), The following discussion reviews several 
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wave-induced current theories and also experimental wave-tank data to determine the 
applicability of such theories. 

This discussion is limited to wave-induced currents only; thus, a number of 
simplifying assumptions are made. Monochromatic, two-dimensional, progressive sur- 
face waves are considered to be moving across the water column. The bottom is postu- 
lated to be flat and impermeable. CorioUs and surface tension forces are neglected, and 
the fluid is considered incompressible, inviscid, and irrotational with a constant, uniform 
surface pressure. A constant mean water-column depth is assumed. Wave-tank tests 
may approximate many of these conditions. 

Given the above considerations, classical wave theory may be used to describe 
theoretically the phenomenon of a wave-induced current. Classical wave theory repre- 
sents a solution of Newton's second law of motion (conservation of momentum) and the 
continuity equation (conservation of mass) subject to boundary conditions prescribed at 
the surface and bottom of the water. Solutions to these equations are found to be 
periodic, taking the form of classical waves. The first-order solution (Airy or linear 
theory) for small ratios of amplitude to wavelength of the waves prescribes a simple, 
sinusoidal waveform in which the fluid particles beneath the surface trace closed orbits 
with the wave passage. Each fluid particle returns to its original position after each 
wave cycle. In infinitely deep water, the orbits are circular; in water of finite depth, 
however, the orbits become ellipses. Higher-order wave theories describe a more com- 
plex situation. 


Stokes' Wave-Induced Mass-Transport Current 

For waves whose amplitudes are not small compared to the wavelength, the equa- 
tion of the free surface is nonlinear. Stokes (ref. 2) found an approximate solution under 
this condition by using a free-surface equation consisting of a series of harmonic func- 
tions. To a second-order approximation, the Stokes solution exhibits the property that 
fluid particles beneath a passing wave trace out orbits which are not closed. After each 
wave cycle, the fluid particles are displaced a distance in the direction of wave advance. 
This nonperiodic drift of the fluid particles (and orbits) has been variously termed the 
mass-transport current, Stokes' drift, and wave-induced current. To a second order, 
Stokes' expression for the mass-transport current profile Ugrp(z) is 


Ust(z) - 


a^ak cosh(j2k(z - h)j 


+ C 


2 sinh'^kh 


( 1 ) 
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where h is the mean water- column depth and C is an arbitrary coefficient. If it is 
assumed that the net mass transport for the water column is zero (assuming constant 
density), 

p \ U(z) dz = 0 (2) 

then 


2 

Q _ -a a sinh 2kh 
4h sinh^kh 

These equations represent the conditions that exist in a closed wave tank where a con- 
stant return flow, given by C, is generated after sufficient time. For deep water 
(kh > 7t), sinh kh may be replaced by e^^/2. When combined, equations (1) and (3) 
become 


^ST,Dp(^) ■ 


a^o-kle 


2kz 



(4) 


Experiments by Mitchum and Mason (ref. 8, p. 32) have confirmed equation (4). Also, 
experiments in reference 6 for 3.2 < kh < 5.5 show reasonable agreement with equa- 
tion (4) in the upper half of the fluid. 

To describe wave-tank experimental current profiles where the wave generator has 
not been running long enough to establish a return flow or to describe idealized ocean 
conditions where a return flow may not exist, the coefficient C is often allowed to be 
zero, thus allowing for a net mass transport. For intermediate waters (kh < tt), this 
assumption yields 


UsT^^^ 


c=o 


a^CTk cosh j2k(z - h^ 
2 sinh^kh 


(5) 


and for deep waters (kh > n) 


^ST,DP^^) 


C=0 


= a^ake"^^^ 


( 6 ) 


for the Stokes mass-transport current profiles. Little data are available to confirm these 
results for finite net mass transport. However, Alofs and Reisbig (ref. 9) have found in 
the movement of lenses of oil on water, under conditions so that 5.50 < kh < 11.0, mean 
surface velocities 35 percent to 150 percent higher than those predicted by Stokes' theory 
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(eq. (6) at the surface (z = 0)). Care was taken to obtain wave -tank data before a return 
flow was set up. ’ 


Stokes' solutions have been obtained to third order (ref. 8, p. 33) and fifth order 
(ref. 10). Other nonlinear theories include cnoidal wave theory (ref. 8, p. 40) and stream 
function theory (ref. 11). These theories represent an effort to match a mathematical 
waveform to the actual one as the wave steepness grows and kh becomes small. All 
these theories, however, neglect the effects of viscosity by assuming an inviscid, irrota- 
tional flow. In intermediate waters (kh < Ji), phenomena occur that are not explained by 
the irrotational theories. Figure 1 shows the experimental mass-transport current pro- 
file for a wave -tank test conducted by Russell and Osorio (ref. 7). Intermediate water 
conditions existed (kh = 1.25), and the data were taken after the wave generator had been 
in operation sufficiently long to establish a return flow and a case of zero net mass trans- 


port for the water column. The wave period T was 1.5 seconds, the wave height H 
was 11.7 cm, and the water- column depth h was 50.8 cm. The Stokes wave -induced 


mass-transport current profile, as computed using equations (1) and (3), is shown as t 
dotted line in figure 1. The theory f'ails^ to predict the observed drift characteristics, 
especially at and near the bottom where a strong positive velocity is present. 


\ 

Longuet-Higgens' Zero Mass-Transport Wave -Induced Current 


Longuet-Higgens (ref; 3) examined Stokes’ assumptions regarding the irrotationality 
of the flow. For a perfect fluid, the flow may have a tangential velocity relative to a 
solid boundary. For a viscous fluid, however, the particles of fluid directly in contact 
with the boundary must possess the boundary velocity, a condition of no slip; at the bottom 
of a wave tank or the ocean this boundary velocity is zero. At a short distance (the 
boundary -layer thickness) above the bottom there exist the oscillatory orbital velocities 
of the fluid particles caused by the wave passage. A strong velocity gradient within the 
bottom boundary exists; hence, strong vorticity. The vorticity may influence the rest of 
the fluid outside the boundary layer by conduction or convection, especially in intermedi- 
ate waters; thus, the flow is rotational. Even if the boundary layer were infinitely thin, 
the vorticity present would be finite. 

In Longuet-Higgens' description, the flow field consisted of a surface boundary 
layer of a thickness 6g and a bottom boundary layer of a thickness 6g in addition to 
an interior (core region) fluid flow (ref. 3). Boundary conditions included zero stress 
at the surface, zero velocity at the bottom, and the imposition of a zero net mass trans- 
port for the entire water column, as shown in sketch (b). 
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Wave propagation 



Longuet-Higgens proceeded with his analysis from an Eulerian point of view. 
Using a perturbation technique and matching the solution for the interior flow to the 
solution at the edges of the boundary layers, Longuet-Higgens showed (ref. 3), for a 
progressive wave where vorticity is assumed to spread only by conduction, that 




a^ak 

4 sinh^kh 




2 cosh 2k(z - h) + A, + A„ + AJ- 


A. = kh sinh 2kh - 3 sinhjl^ _ 3 
1 2 kh 2 


A2 = -4kh sinh 2kh 


> ■ ( 7 ) 


Ag = 3kh sinh 2kh 


3 sinh 2kh ^ 9 
2 kh 2 




This equation applies to the interior of the fluid. , Separate equations (ref. 3) exist for,, 
the mass -transport current profiles within the thin boundary layers. ^ However, equa- 
tions (7) describe the core region drift flow usually measured in wave-tank experiments. 
Shown in figure 1 is the Longuet-Higgens solution for the, mass-transport current 
profile based on the previously described wave conditions. It is. evident that -a better 
qualitative fit to the experimental data curve is obtained. This observation is, especially 
true at and near the bottom (outside the bottom boundary layer) where Stokes predicts a 
negative velocity, whereas the Longuet-Higgens theory predicts the strong forward flow 
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actually observed. It should be noted that within the bottom boimdary layer (in this case 
0.7 mm thick) the actual velocity would be zero at the bottom wall as is predicted in the 
bottom boimdary-layer solution of Longuet-Higgens (ref. 3). Stokes' theory, which 
assumes inviscid flow, predicts a slip flow. 

The conduction solution of Longuet-Higgens (eqs. (7)) is strictly valid only for wave 
amplitudes very much smaller than the boimdary-layer thicknesses (which themselves 
are small), that is, a^ 6g - a/6g « 1. JHowever, wave-tank tests by Russell and Osorio 
(ref. 7) show the equation to be of use for waves of considerable amplitude. For the 
particular example of figure 1, a is 0.058 meter and Sg 6g =: 0.7 mm yielding 
a/6g ~ = 83. Generally, Longuet-Higgens' theory, based on the experimental evi- 

dence of Russell and Osorio (ref. 7) where 
favorably wave-induced drift profiles in waters where 0.7 < kh < 1.5. For kh < 0.7 
the waves begin to exhibit an instability invalidating the Longuet-Higgens perturbation 
assumption. In deeper waters (kh > 1.5), the assumed surface condition leads the 
Longuet-Higgens solution to an unbounded condition near the surface (ref. 12). 

Longuet-Higgens' analysis describes the wave-induced mass -transport current 
profile for a viscous, intermediate water condition of zero net mass transport of the 
water column. Closed wave-tank experiments have confirmed essential features of this 
theory for zero net mass transport. However, for cases of a nonzero mass transport, 
as in a long wave tank where the wave generator has not been running long enough to 
establish a return flow, or in the ocean, it is unlikely that a zero net mass transport 
caused by waves exists. Therefore, the use of a zero mass transport, Stokes' equation 
(eq. (1) or eq. (4)) or Longuet-Higgens' equations (eqs. (7)), would not be valid. Instead, 
Stokes' nonzero mass-transport theory (in deeper waters where viscous effects may not 
be important for most of the flow profile, eq. (5) or (6)) or an intermediate- and shallow- 
water viscous theory with nonzero mass transport should be used. To be of most value, 
any such theory would require confirmation by experiment. 


V^S ~ V^B ^ 0(100), appears to predict 


Unliiata and Mei's Nonzero Mass- Transport Theory 

In a paper by Unluata and Mei (ref. 4), boundary -layer arguments were incorpo- 
rated into the Lagrangian form of the equations of motion where the free surface may be 
treated as a fixed and known quantity (unlike the Eulerian Longuet-Higgens theory where 
complicated coordinate transformations of the free surface are required). Again, the 
fluid is divided into surface and bottom boundary layers with a core (interior) region. By 
perturbation analysis the second-order terms yield the following solution for the wave- 
induced current 
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( 8 ) 





(z2 - h2) + 

- + (l - -'j kh sinh 2kh + ^ cosh 2k(z - h) 

' sinh'^kh 

[J^ \ h) 2 J 


and for the wave-induced net volume transport 

V = + — sinh 2kh + — sinh 2kh\ (9) 

^ 8p sinh^kh v ^ / 

In these equations, is a time-averaged, constant pressure gradient associ- 

ated with the wave-induced current in the coordinate system of the fluid particles. 
According to Unluata and Mei, two limiting cases are possible; (1) when = 0, the 
case of zero net volume transport induced by waves, equations (8) and (9) confirm the 
Eulerian approach of Longuet-Higgens (eqs. (7)); (2) in an infinite wave tank or idealized 
ocean where steady state has been reached = 0, and a maximum net volume trans- 

port of the water column caused by waves occurs. Also, according to Unluata and Mei, 
the equations are most appropriate for 0.7 < kh < 1.5 when compared to the Russell and 
Osorio experimental results. No discussion of possible intermediate cases of net volume 
transport is given in reference 4. 


Experimental Data 

Few experiments have been performed to study wave-induced drift profiles under 
the conditions of a nonzero volume transport. The experiments of Alofs and Reisbig 
(ref. 9) were conducted in a closed wave tank under deep water conditions (kh > 5.5) and 
controlled conditions of a nonzero volume transport. Only surface drift velocities 
caused by waves were measured, but these velocities proved greater by 35 percent to 
150 percent than predicted by the Stokes nonzero mass-transport theory. Oil lenses 
were used as surface drifters which may move at a different velocity than the water 
beneath them and influence the surface boundary- layer conditions also. Thus, it is diffi- 
cult to assess these results when they are compared to a clean water condition. Some 
experiments, reported in reference 6, were conducted in intermediate waters (kh < tt) 
before a return flow was set up, thus producing conditions of a possible nonzero volume 
transport caused by waves; but the results show a considerable scatter in the data. It is 
suspected, as reported in reference 6, that the scatter may be a result of the failure to 
achieve a stable waveform at the data station when the data were taken. Since the wave- 
induced mass-transport current is a secondary effect, the error introduced by an unstable 
waveform may be considerable. Thus, little usable information is available in cases of a 
nonzero volume transport caused by waves. 
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In wave-tank experiments, the usual measured data include (1) the wave and tank 
properties, and (2) the surface or bottom drift velocities, or the entire observed wave - 
induced drift profile. Small floats, neutrally buoyant particles, or dye streaks are the 
usual method of tracking water movements. The Unliiata and Mei theory (eqs. (8) . 
and (9)) suggests that the second-order, time-averaged pressure gradient P 2 ^ is a 
necessary property for a comparison of the predicted and measured drift profiles. In 
practice, may be a difficult parameter to measure with sufficient accuracy. It is 

suggested, therefore, that an alternate form of the equation for the wave -induced mass- 
transport current profile may be more suitable for experimental purposes. A reexami- 
nation of the Longuet-Higgens analysis, explained in the next section, shows that such an 
equation may be determined. 


LONGUET-HIGGENS' NONZERO MASS -TRANSPORT EQUATION 
FOR EXPERIMENTAL APPLICATIONS 


Longuet-Higgens, in his paper on mass transport caused by waves (ref. 3), alluded 
to the possibility of a nonzero mass transport for the water column as a whole, and even 
suggested the form any additional terms, added to equations (7), might exhibit. However, 
he did not explicitly derive and state such an equation. An equation of this form would 
help describe those experiments conducted in wave tanks under nonzero mass- transport 
conditions. Such an equation would require only those wave parameters usually measured 
in these experiments including wave height, period, wavelength, the wave-induced current 
profile, and water depth. . 

To develop such an equation, Longuet-Higgens' original derivation was followed to 
the point where he assumed a case of zero mass transport (or zero volume transport for 
a constant density fluid) for the entire water column. At that point, an arbitrary nonzero 
total volume transport was assumed to exist. The derivation of the final form of 
the theoretical equation is presented in the appendix. Equations (7) are modified to 




^LH,MT^^) ■ ^LH^^) ^ 2 "K 


1 - 


( 10 ) 


for cases of nonzero volume transport of the water column caused by waves. It is also 
noted, from material presented in the appendix, that 


V. 


X 


— 

Jo 


( 11 ) 
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The previously stated assumptions regarding the Longuet-Higgens theory are assumed 
to apply, with equation (10) most applicable in the range 0.7 < kh < 1.5. 

Equations (10) and (11) represent an arbitrary perturbation on the Longuet-Higgens 
zero mass-transport profile to account for a nonzero mass transport. Experimental 
data, as yet unavailable, must be used to validate this new approach. 

A number of general design criteria exist in the planning of an experiment in a 
wave tank to study conditions of nonzero mass- transport wave-drift profiles. Wave 
amplitudes must be kept small so that the quantities a/Sg a/6g are small although 
the Longuet-Higgens theory was shown by Russell and Osorio to be useful even when 
ay^6g = a^6g = 0(100). Surface contamination of the water must be avoided. Huang 
(ref. 12) has shown such contamination to have a pronounced effect on the surface bound- 
ary condition and the resultant wave-drift profiles. Alofs and Reisbig (ref. 9) discuss a 
method whereby such effects of surface contamination are minimized. The method 
involves using a fan at the generating end of the wave tank to move the surface contami- 
nants to the other end of the tank where they may be trapped or removed. 

The time required for a backflow current to develop after wave generation has 
started in a wave tank is of the order of £/U(0) where £ is the length of the wave 
tank and U(0) is the surface wave-induced drift velocity (ref. 4). Longuet-Higgens 
(ref. 3) indicates that the time required to establish a steady state with regard to the 
conduction of vorticity is of the order of h^/ 1 >. To produce a condition of maximum 
wave-induced volume transport mAX’ the wave generator should have been running 
long enough to establish a stable waveform and a steady state with regard to the conduc- 
tion of vorticity. The wave generator should not, however, run so long that a return 
flow can be felt at the data measuring stations. The above discussion indicates that max- 
imizing the wave-tank length £ and minimizing the wave-tank water depth h would 
produce the desired effects. Then, as time passes and a return flow develops, the net 
volume transport would decrease from its maximum value MAX '^^itil, for a fully 
developed backflow, the net volume transport is zero (V^ = 0). At the far end of the 
wave tank, a sloping beach with wave -absorbing devices and materials should be used to 
prevent wave reflections and standing waves which would not satisfy the purely progres- 
sive wave assumption made in this analysis. 


Under the experimental procedures outlined, an experimenter might obtain wave- 
induced current profiles U(z) (for a particular set of wave conditions) at successive 
times corresponding to varying backflow current conditions. Just after wave generation 
has begun (zero backflow current) a maximum net volume transport MAX 
with the corresponding wave-induced drift profile Ug^P time, a backflow 

current develops, and the wave -induced current profile proceeds through a series of 
stages corresponding to successively smaller net volume transports Vj^. Finally, at a 
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sufficiently long time after the beginning of wave generation, a fully developed backflow 
current exists corresponding to zero net volume transport (Vx = 0) and to the waver 
induced current profile Uj,^p(z). Equations (10) and (11) may be rewritten as 


and 


jT t rw\ TT ^ ]\^ AX 



^x,MAX Jq ’^EXP,MT^^^ 


( 12 ) 


(13) 


where the experimental values of wave-induced current profiles at maximum net volume 
transport Ug^P zero net volume transport UgxP^^^ used in the calcu- 

lations. Confirmation of the perturbation on the Longuet-Higgens zero mass -transport 
profile would exist if 


" ^EXP,MT^^^ 


(14) 


is obtained for the case of maximum net volume transport. Additional checks may be 
obtained by using the intermediate stage wave-induced current profiles corresponding to 
net volume transports 0 ^ = Vj^. equations (12) and (13). 

For equation (10) to be a predictive technique, however, the net volume trans- 
port Vjj must be parameterized. This process would, of necessity, involve a series of 
wave -tank tests to find the dependence of on the wave and wave-tank test conditions. 

Thus, in general. 


= V^(a, L, T, h, .) 


(15) 


By dimensional-analysis and data-analysis techniques it should be possible to reduce 
the number of independent variables and to find a relation for based on the wave- 

tank test data. The theoretical relation for MAX Unluata and Mei (ref. 4) for the 
limited range of applicability (0.7„< kh < 1.5) would be helpful in this analysis. 

In the absence of experimental data obtained under conditions of a nonzero volume 
transport caused by waves, it is informative first to analyze parametrically some typical 
wave-tank and ocean-swell cases to examine: (1) the variation in wave-induced current 
profiles between extremes of zero net volume transport and maximum net volume trans- 
port and (2) typical magnitude ranges for the wave-induced currents. In this manner, the 
results of the next section may shed some light on what might be expected in an experi- 
mental program. 
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RESULTS AND DISCUSSION 


The example of figure 1 has been examined for a range of possible volume trans- 
port values. Again, the particular wave conditions for this experiment were: 


Water depth, h, cm . 
Wave period, T, sec 
Wavelength, L, cm . 
Wave height, H, cm 
kh . . . 

^As ~ V^B 


50.8 

1.5 

256 

11.7 

1,25 

83 


The experimental data presented in figure 1 represent a condition of zero net volume 
transport. Russell and Osorio (ref. 7) took this type of wave-tank data after the wave 
generator had been in operation for hours, when necessary, to achieve steady backflow 
conditions. The Unliiata and Mei equation (9) was used to obtain a maximum theoretical 
volume transport. According to Unliiata and Mei, this maximum occurs when ^2a. ~ ^ 
corresponding to an infinite wave tank. Equation (9) can be written for this case as 


^x,MAX 


^ ^ 2kh + 
sinh^kh\4 2 


1 

4kh 


sinh 2kh 


( 16 ) 


The range of possible volume-transport values chosen for this parametric analysis 
included: (1) a solid wall wave -tank end (Vx = 0), (2) various degrees of wave-tank end 
"porosity" (Vx = 0.25 V^^mAX> ^x = 0. 50V^^jVlAX’ (3) an 

open-ended, infinite length channel (Vx^mAx)- With these values and the previously 
stated wave conditions, the Longuet-Higgens nonzero mass -transport equation (10) 

calculated from eqs. (7)j. The results are plotted in figure 2. 

In figure 2 the case of Vx = 0 corresponds to the Longuet-Higgens zero volume - 
(or mass-) transport case originally presented in figure 1. However, as more net volume 
transport caused by waves is assumed (i.e., a smaller backflow), the effect is to shift the 
entire velocity profile forward to more positive values, thus rotating the curve about the 
point of interior flow near-bottom velocity. This represents a parabolic velocity distri- 
bution added to the original Longuet-Higgens equations (eqs. (7)) that vanishes at the 
bottom boimdary layer and has a zero velocity gradient at the free-surface boundary layer 
as Longuet-Higgens had suggested (ref. 3). Finally, there is the extreme theoretical 
case of the infinite or open-ended channel where the volume transport is a maximum. A 
substantial wave-induced current profile, everywhere positive in the direction of wave 
advance, is present. 


was used 


fwith 


Ulh(2) 
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A case was also examined of a typical ocean swell, originally considered in refer- 
ence 1. The particular wave conditions for the swell were: 


Water depth, h, m 18.3 

Swell period, T, sec 11.0 

Swell wavelength, L, m 132.31 

Swell height, H, m . . . 1.83 

kh 0.87 


The value of kh in this example (0.87) is within the intermediate water limits best 
applicable for the use of the Longuet-Higgens nonzero mass -transport equation (10) 

(0.7 < kh < 1.5). Examination of the i-atio a/6 , is important to check the validity of 
using the Longuet-Higgens theory. To be strictly valid, a/6g ~ a/6g should be much 
less than one, but Russell and Osorio showed the theory to predict reasonable results 


when ay6g 
defined as 


order of 100. The boundary -layer thickness 6 


IS 


5 = 



(17) 


For water the kinematic viscosity u is about 0.01 cm^/sec and for this swell wave 
a = 2tt/T = 0.571 sec"l. This yields a boundary -layer thickness 6g ~ 6g of 0.187 cm. 
The ratios ay^6g and ay'Og are calculated to be 489. These values are almost six 
times the value calculated for the wave-tank tests of Russell and Osorio. The question 
posed is whether or not the Longuet-Higgens theory may be applicable for these large 
values of a/6. For this reason, extreme caution must be exercised in the interpreta- 
tion of the results of any calculations made for this swell condition. Therefore, until 
experimental data are in hand, the following discussion is largely conjectual. 


The swell height chosen (1.83 m) for this case is a winter statistical value for the 
Chesapeake Bight region so that 55 percent of the time the swell was at or below this 
value and 45 percent of the time above this value. Thus, the swell height may be con- 
sidered typical. The Stokes mass -transport current calculated from equation (5) is shown 
as the dashed line in figure 3. The solid curves represent Longuet-Higgens' mass- 
transport current profiles using equation (10), with volume-transport values of = 0, 

Vx = 0.25 Vx,MAX’ Vx = 0.75Vx,MAX’ ^x = Vx^mAX where Vx^mAX was again 
determined from Unliiata and Mei's theory (eq. (16)). As figure 3 shows, a considerable 
range of current profiles is possible depending upon the volume -transport conditions. In 
all cases, the velocity near the bottom predicted by the nonzero volume-transport Longuet- 
Higgens theory is greater than that determined from Stokes' theory which neglects vis- 
cosity. For zero volume transport, a negative midlayer drift exists; whereas, for most 
positive values of volume transport Vx the drift profiles are everywhere positive. The 
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surface drift increases with increasing volume transport. In the cases of small net vol- 
ume transport, the wave-induced surface current is on the order of 1 to 2 cm/sec, an 
amount considerably less than the static circulation of 5 to 20 cm/sec generally found 
along the Atlantic continental shelf (ref. 1, p. 24). At the higher values of net volume 
transport, however, the wave-induced surface current is of the order of 10 cm/sec or 
the same order as the static circulation on the continental shelf. 

If these predictions are correct, the surface and midlayer drifts of pollutants, 
sediments, and buoys may be affected by the wave-induced current. In the case of buoys, 
large drogue plates or parachutes are set for a predetermined depth for the purpose of 
measuring the current at that depth. When the drogue plate or parachute drag area is 
much greater than that of the surface float and interconnecting cables, then the buoy 
drifts essentially with the current at the drogue plate or parachute depth. The surface 
float marks the position for quick determination of drift. Examination of figure 3 shows, 
for example, that for a particular depth on the continental shelf where the drogue is at a 
depth of 0!3 of the water-column depth, for zero volume transport the drift caused only 
by the swell-induced current is opposite to the direction of swell propagation (-1 cm/sec); 
however, for a case of maximum volume transport, the buoy drifts quickly in the direc- 
tion of swell propagation (8 cm/sec). Thus, in the presence only of swell-induced cur- 
rents, by this theory completely opposite buoy drifts are predicted dependent upon the 
net volume transport that exists. Further, in the interpretation of a known buoy drift 
where it is desired to subtract out, in some analytical manner, the influences of mass- 
transport currents (as well as tidal currents, wind stress currents, etc.) to obtain a 
residual "mean current" (ref. 8, p. 336), it is important to determine properly the total 
volume transport caused by waves to obtain the best estimate of the wave-induced mass- 
transport current profile. 

Detailed, open-water continental shelf experiments which can be used to test the 
wave-induced mass-transport current calculations using equations (10) and (11) are not 
yet available. Wave -tank tests, however, appear to be a reasonable alternative for 
testing the Longuet-Higgens nonzero mass -transport theory based on the design criteria 
outlined previously. These tests would also aid in parameterizing the net volume trans- 
port so that equation (10) may be used eventually as a predictive equation for future 

studies. 


CONCLUDING REMARKS 

4 

This paper presents a review of some wave-induced mass-transport current 
theories both with and without a net mass (or volume) transport of the water column. 
A theoretical relationship based on the Longuet-Higgens theory was derived for the 
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wave-induced mass-transport current profiles in intermediate waters with viscous flow 
for cases of nonzero volume transport. The theoretical relationship is presented in a 
form usefiil for experimental applications. Some design criteria for experimental wave- 
tank tests using this relation are also presented. In the absence of such experimental 
data, sample parametric cases for monochromatic waves for a typical wave-tank test and 
a typical ocean swell were assessed using the derived relation. A maximum wave- 
induced volume transport was determined from an equation of Unliiata and Mei as an 
upper limit for the parametric cases. 

Calculations indicate that substantial changes in the wave-induced mass -transport 
current profiles exist dependent upon the assumed net volume transport. The maximum 
volume transport, corresponding to an infinite channel or idealized ocean condition, pro- 
duces the largest wave-induced mass-transport currents. The surface and midlayer 
water-column currents caused by waves increase with increasing net volume transport. 
These results suggest that wave-induced mass-transport currents may have a consid- 
erable effect upon pollution and suspended sediment transport and buoy drift. Also, the 
parametric analysis indicates that a more detailed examination of wave-induced mass- 
transport currents based on extensive wave-tank testing is necessary to substantiate the 
Longuet-Higgens nonzero mass-transport theory and to parameterize the wave-induced 
net volume transport. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
November 3, 1975 
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APPENDIX 


AN EQUATION FOR WAVE -INDUCED MASS -TRANSPORT CURRENT 
PROFILES FOR INTERMEDIATE WATER, VISCOUS FLOW 
WITH NONZERO VOLUME TRANSPORT 

Longuet-Higgens (ref. 3) examined the problem of wave -induced currents for a 
fluid of small viscosity. His analysis included solving for the stream function of the 
interior flow (between the edges of the surface and bottom boundary layers), but the 
analysis used the conditions at the interface of these layers as necessary boundary con- 
ditions. By assuming a zero net mass transport of the water column, Longuet-Higgens 
was able to obtain a single solution for the stream function. Longuet-Higgens alluded to 
the possibility of a nonzero net mass transport caused by waves and suggested the form 
any additional terms added to the stream function would exhibit. This appendix is an 
explicit statement of his inferences. An equation is derived for the wave-induced veloc- 
ity profile in a form that would be useful for experimenters studying wave-induced cur- 
rents. It is shown that the term added to the original Longuet-Higgens zero mass- 
transport equation is, in fact, of the form suggested by Longuet-Higgens. 

This appendix closely follows the Longuet-Higgens derivation for a zero net mass 
transport, and the reader may consult reference 3 for additional details. Only mono- 
chromatic, two-dimensional, progressive waves are assumed to be moving through the 
water column in a manner similar to waves in a wave tank or swell in the ocean. A con- 
stant density fluid is assumed. The net mass transport, then, is simply the product of 
the constant density and the net volume transport. The boundary layers at the surface 
and bottom of the fluid 6g and 6g are very small compared to the depth of the interior 
flow. Thus the distances from the surface to the top and bottom interfaces may be 
written 


z = 6g = 0 
z = h - 6g = h 

In this instance equations (291), (286), (287), (288), and (289) of reference 3 may be 
rewritten,, using the notation of this paper, as 

vV(z) = v4 sinh 2k(z - h) (^2) 

. . 4 sinh^kh 


(Top interface) 


(Bottom interface) 


(Al) 
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e 



5oka^ 

4 sinh^kh 


(A3) 


^2/ j _ -2ok^a^ sinh 2kh 

(A4) 

sinh^h 

= 0 

(A5) 

~ Constant 

(A6) 


Equation (A2) relates the stream function for the interior of the fluid to the wave and 
mean water-column depth conditions and is a function of z. Equation (A3) is, a boundary 
condition at the edge of the bottom boundary layer. The equation represents, in fact, the 
horizontal wave -induced drift velocity at that point. The quantity e is a small dimen- 
sionless quantity of the order of a/L. Equation (A4) represents the gradient of the 
wave-induced drift velocity at the surface boundary layer. The stream function \}/ is 
constant at the boundaries. At the surface, this constant may be chosen to make \p 
vanish, yielding equation (A5). The stream function at the bottom boundary layer is also 
a constant. Longuet-Higgens examines this and, for the case of zero net wave-induced 
mass transport, sets equation (A6) equal to zero. As Longuet-Higgens indicates, however, 
this condition may be relaxed to allow a nonzero net mass transport. As e^i//(z) is the 
stream function for the mass transport velocity MT^^^ (ref. 3), then 

Iz " ^LH,MT^^) ' 

Integrating equation (A7) over the water column and considering equation (A5) results in 

9 

The value of e (;//) is the wave -induced mass -transport velocity profile integrated 
over the water column. This value is equivalent to the net volume transport (per unit 
crest width of the wave). If an arbitrary value for the integral (e.g., Vx) is prescribed 
instead, then equation (A6) becomes 

e^(iL>) , = V 

^z=h X 

where is the volume transport per unit length of the wave crest. 
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To solve equations (A2) to (A5) and (A9), a solution of the form 


— [sinh 2k(z - h) + R(z)] 
4 sinh^kh 


-(AlO) 


is assumed as proposed by Longuet-Higgens for the case of zero net volume (mass) trans- 
port. It is then noted that 


aa^k 


LH,MT az 4 sinh^kh 


2 cosh 2k(z - h) + i 

k 9z 


(All) 


where Uljj indicates the expanded Longuet-Higgens equation to account for a 

nonzero volume (or mass) transport. 

Equation (AlO) may be differentiated as required, and the boundary conditions at 
z = 0 or z = h substituted as needed, using equations (A2) to (A 5) and equation (A9). 
After the necessary algebraic operations, 

(A12) 
(A13) 

(A14) 

[r(z)"]^_q = sinh 2kh (A15) 

and 


d R(z) ^ Q 
dz^ 


dR(z) | 

dz 


= 3k 


lz=h 


\ d\z) \ 

dz‘^ 


= -4k sinh 2kh 


lz=0 


[r(z)1 , (A16) 

^ -*z=:h a^CT ^ 

are obtained. Again, equations (A12) to (A15) are the same as presented by Longuet- 
Higgens, but equation (A16) is different from zero, thus the effects of a net volume trans- 
port are indicated. 

From equation (A1 2), R(z) must have the form 

R(z) = “ 3 ^ + *^4 (A17) 
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where a-^, 0 . 2 , and are functions only of the wave and water-column depth 

conditions and the assumed value of volume transport. Using equation (A17) in equa- 
tions- (A13) to (A16), four equations result in four unknowns: a-^, a 2 , ^3, and a^. 
Solution of these equations yields 


Q- = -2 V +-^ + — sinh 2kh -f sinh 2kh 

^ h^a^a ^ 2h2 h 2h3 


Ofg = -2k sinh 2kh 


= + lk\ - I-) sinh 2kh + 6 si^^kh „ 

2 2 \ 2hJ a2ah 


0!^ = sinh 2kh 


(A18) 


Differentiating equation (A17) with respect to depth yields 


+ 2a,^z + o!„ 
dz 1 z 6 


(A19) 


Equations (AH), (A18), and (A19) may all be combined to obtain 


where 


U 


a^ak 


(7\ - 


kh 


2 cosh 2k(z -h)+B-^ + Bgl^j -1- Bgl J 


(A20a) 


Bi = kh sinh 2kh - | - I + ^ ^ 


2 kh 


a'^akh 


B2 = -4kh sinh 2kh 

B„= 3kh sinh 2kh + ^ 2kh ^ 9 _ 6 sinh^kh y 

o 9 ‘ ' 


2 kh 


a'^CTkh 


j 


(A20b) 


The basic difference between this expanded form of the Longuet-Higgens current profile 
equation for nonzero volume (or mass) transport and the original Longuet-Higgens equa- 

6 sinh2kh 

tion (ref. 3) is the additional term s V„ which appears in both the B^ and 

a^akh ^ ^ 
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coefficients. This term represents changes as a result of a net volume transport 
caused by the waves. Finally, equations (A20a) and (A20b) may be rewritten as 






1 - I- 


(A21) 


where Uljj(z) is the Longuet-Higgens wave-induced current profile equation for zero 
net volume transport (ref. 3). Longuet-Higgens indicated that a term representing a 
parabolic velocity distribution, vanishing at the bottom and with zero velocity gradient 
at the free-surface boundary layer, would be added to the current profile equation for 

Ulh(z). 


The term containing in equation (A21) is of this form. The net volume 


transport is an unknown quantity which must be parameterized, if not theoretically. 


then by experiments, to obtain unique solutions of equation (A21). 
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Wave-induced mass transport current, cm/sec 
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Figure 1.- Wave -induced mass -transport current plotted against nondimensional depth for zero 
net volume transport = oj . Surface and bottom boundary layers are not shown. 




Longnaet-Higgens' mass transport current, cm/sec 



Figure 2.- Longuet-Higgens' mass- transport current plotted against nondimensional depth 
^with and without net volume transport of water column^ using equation (10). Surface 
and bottom boundary layers are not shown. 









Swell-induced mass transport current, cm/sec 
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Surface and bottom boundary layers are not shown. 
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